Coherent Synchrotron Radiation for Laminar Flows 
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We investigate the effect of shear in the flow of charged particle equilibria that are unstable to the 
Coherent Synchrotron Radiation (CSR) instability. Shear may act to quench this instability because 
it acts to limit the size of the region with a fixed phase relation between emitters. The results are 
important for the understanding of astrophysical sources of coherent radiation where shear in the 
flow is likely. 



I. INTRODUCTION 

When the wavelength of synchrotron radiation emit- 
ted by a bunch of relativistic particles is comparable 
to the size of the bunch the particles may radiate co- 
herently. The coherent emission produces significantly 
greater power than the incoherent emission. Coherent 
synchrotron emission is seen from different astrophysical 
source but most notably from the rotationally powered 
radio pulsars (e.g., [l|). In particle accelerators the coher- 
ent synchrotron emission is usually undesirable because 
it causes very rapid energy loss. In this context we would 
also like to mention efforts to study coherent curvature 
radiation as a pulsar emission mechanism in the labora- 
tory 0. 

Several effects may act to stabilize a system unstable 
to the CSR instability. In 0] the authors analyzed the 
influence of a small energy spread in a beam of charged 
particles in approximately circular motion. A distribu- 
tion function with a single value of the canonical angular 
momentum was considered. The radial width of the beam 
was given by the amplitude of the betatron oscillations 
which is non-zero for a non-zero energy spread. It was 
shown that the decohcrence introduced by the betatron 
oscillations leads to a characteristic frequency spectrum, 
whereas the dependence on the Lorentz factor and the 
number density remains unaffected. 

In this paper we relax the assumption of zero spread in 
the canonical angular momentum P^ of the equilibrium 
distribution. In general this leads to shear in that the av- 
erage angular velocity becomes dependent on the radius 
of the orbit. In practice the requirement of a small spread 
in the canonical angular momentum may be harder to 
satisfy than a small energy spread. Shear is expected in 
astrophysical sources of coherent emission 

Shear itself can be the cause of instabilities, for exam- 
ple, the "diocotron instability" [fj, but this is not the 
focus of the present paper. In the case of CSR it is 
reasonable to expect that the shear acts to stabilize the 



CSR instability due to particles with different radii "slip- 
ping away" . Using a linear perturbation analysis of the 
fluid equations for a laminar Brillioun flow we show that 
even with a spread of the canonical angular momentum 
AP^ > 0, previous results for equilibria with constant 
P$ can be recovered treating the plasma as a relativis- 
tic cold fluid. Computer simulations of CSR emission in 
Brillioun flows p are compatible with this picture. With 
the cold fluid approximation adopted here, the stability 
depends on the number density and the angular velocity 
which depends on the radius. The azimuthal velocity is 
approximately equal to the speed of light. 

The problem of CSR emission has been investigated by 
several authors. Goldreich and Keeley Q considered the 
stability of a charge distribution whose motion is con- 
fined to a thin ring with the particle motion being one- 
dimensional. This calculation led to confusion as to how 
the proposed CSR mechanism works in detail (cf. [8|). 

As proposed in 0] the CSR instability is related to 
the classical negative mass instability [1, [l(| [LH G3 in 
the sense that an increase in particle energy leads to a 
decrease in its angular velocity. While the classical nega- 
tive mass instability is caused by the Coulomb part of the 
electromagnetic potential, the CSR instability is caused 
by the radiation field. The negative mass effect is not im- 
mediately apparent in a one-dimensional treatment based 
on conservation of energy and charge. Heifets and Stu- 
pakov [l]| effectively built in the negative mass effect by 
hand having a constrained radius and particle energy. 

The treatments [3, 0, El give the growth rate in the 
absence of an energy spread. However, considering a 
non-zero energy spread requires a truly two-dimensional 
model. Larroche and Pcllat investigated the effect of 
steep boundaries in the particle distribution function [l4| . 

Section II describes the assumed models, and Section 
III approximate solutions of the equations. Section IV 
derives a dispersion relation for the considered perturba- 
tion. Section V discusses the results. 



II. THEORY 
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We consider two cases illustrated in Figure 1. The 
first case (a) corresponds to a thin layer of relativistic 
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FIG. 1: Panel (a) shows the geometry of a relativistic E- 
layer for the case of a uniform external axial magnetic field. 
Panel (b) shows the second case of a relativistic E-layer in an 
external toroidal magnetic field and an external radial electric 
field. 



We consider weak layers in the sense that the 'field 
reversal' parameter 

C = - — / drp e v (3) 

is small compared with unity, C 2 1. Under this condi- 
tion equation [1] gives fl — —qB e /(m e -f). Here, we have 
assumed that the layer exists between r\ and r^. We 
also consider that the Lorentz factor is appreciably larger 
than unity in the sense that -f 2 1. Furthermore, we 
consider radially thin layers 

Ar 

S = — < 1 . (4) 

We consider general electromagnetic perturbations of 
the electron layer with the perturbations proportional to 

f a (r) exp(iTO(/> — iut) , (5) 

where a — 1,2,.. for the different scalar quantities, m = 
integer, and to the angular frequency of the perturbation. 
Thus the perturbations give rise to field components SE r , 
SE<p, and SB Z . The perturbed equation of motion is 



electrons gyrating in a uniform external magnetic field, 
an "E-layer." The second case (b) corresponds to the 
electrons moving almost parallel to a toroidal magnetic 
field. 



A. Case a 

We consider a laminar Brillioun type equilibrium of 
a long, non-neutral, cylindrical relativistic electron (or 
positron) layer in a uniform external magnetic field B e = 
B e z, where we use a non-rotating cylindrical (r, (j>, z) co- 
ordinate system. The electron velocity is v = t>0(r)0 = 
rrt{r)4> The self-magnetic field is in the z— direction 
while the self-electric field is in the r— direction. The 
radial force balance of the equilibrium is 



vB) 



(1) 



where 7 = (1 — w 2 ,) -1 ^ 2 is the Lorentz factor with veloci- 
ties measured in units of the speed of light, B = B e + B s 
is the total (self plus external) axial magnetic field, E is 
the total (= self) radial electric field, and q and m e are 
the particle charge and rest mass. We have 

IdirE) t dB 

= 47rp e , — = -4:irp e v , 2 

r ar ar 

where p e (r) is the charge density of the electron layer. 



— (c5E + v x SB + Sv x B) , 



(6) 



where the deltas indicate perturbation quantities. This 
equation can be simplified to give 



—ijAu) 
jil + ( 7 fir)' + £-B 



^n(i + ^)-^-B 



Sv r 
8v<i, 



q 



6E r + vSB z 
5E 6 



(7) 



where the prime denotes a derivative with respect to r, 
and 

Aw(r) = lu — mQ(r) 

is the Doppler shifted frequency seen by a particle rotat- 
ing at Q. We also define the dimensionless quantity 



Aw 



Au 



(8) 



which will turn out to be useful later. 

Using the equilibrium equation [T] and the condition 
( 2 <C 1, the matrix in equation [7] is approximately 



V 



—i^Aiv — 7 3 £1 
7 3 (fir)' -ij 3 Au 



(9) 
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We have used the fact that ftfir)' = 7 3 (Or)'. For C 2 < 1 
we have (Or)' = O/7 2 and il' = — v 2 £l/r. In the absence 
of shear the latter quantity would be zero. Consequently 



det(Z>) = 7 4 (ft 2 - Aw 2 ) 
Inverting equation [7] gives 



Sv r 
ami 



97 



2 det(X>) 



97 

m e det(2?) 



iAujSE,/, - Q(SE r + v5B z 



(10) 



, (11) 



(12) 



Here, qp e j 3 £l/ (m e det(X>)) has the role of the distribution 
function of angular momentum [l5| . 



B. Case b 

Here we consider an equilibrium with the same number 
density and velocity profile as before but with different 
external fields as shown in Figure [TJd. Instead of an ex- 
ternal magnetic field in the z— direction we consider an 
equilibrium with an azimuthal magnetic field acting as 
a guiding field and a radial electric field. The latter is 
included in the equilibrium condition and therefore does 
not enter the linearized Euler equation. would only 
enter if we considered motion in the axial direction and 
non-zero axial wavenumbers. 

Thus, we obtain the matrix V again without the Bq 
terms, i.e. for 7 S> 1 



V 



with 



We obtain 



-i^Alu — 7 3 f2 
27O -i7 3 Aw 

det(D) = 7 4 (20 2 - Aw 2 ) 



Sv r 



and 
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m e det(X>) 



97 
e det(X>) 



iAuj(SE r + v5B z ) + £16E 4 



iAuSE^ - 2Q{5E r + vSB z ) 



(13) 



(14) 



(15) 



(16) 



Such a configuration is a more realistic possibility for the 
magnetosphere of a radio pulsar. The conclusions for the 
two configurations do not differ significantly, and we will 
proceed analyzing case a. 



III. APPROXIMATE SOLUTION 

Using equations (fTTj) and (TT2"]) the linearized continuity 
equation gives 



i Aujdp 



r 

1 d 



qipo 



m e det(X>) 



{-iAujSE^ - ft (SE r + v^SBz)) 



97 3 Po 



r dr [ m e det(2?) 



(-iAw (SE r + 



) + nsE^ 



(17) 



where we used SJ r = poSv r and SJ^ — v^Sp + poSv^,. 
We consider conditions where SE r + v l f,8B z terms can be 
neglected. The sufficient conditions are 



Alj 



n 



< 1 



and 



\SE r + vSB z \ < 



Alu 



1^1 



(18) 



(19) 



We estimate the relative magnitude of the field compo- 
nents using the Maxwell equations. From Faraday's law 
we obtain the relation 



D r (8E$) — ik(j,5E r = iu;5B z 



(20) 



where D r (...) = r~ 1 d/dr(r...) « ik r assuming the radial 
dependence exp(ifc r r) for the perturbed quantities as well 
as thin layers with 5 <C 1. We have 

krSEtj, = kff, (SE r + v<p5B z ) + kj,v$Au)8B z (21) 

v^SQ and 5E<f, = 



With 5B Z = DrSAj,, SA^ 



5E r + v^SBz 



Inequality (|19p turns into 

r Aw 
k r <C m 



m Alj 



6E&2) 



(23) 



where k r = k r ro. We will always assume radial wavenum- 
bers k r that are sufficiently small such that the latter 
condition is met. A consequence of inequality (fTS]) is 
that 



Au) 



(24) 



Making use of approximation 1191 equation [TT] can be 
entirely written in terms of 6Em. Note that the resonant 
term due to det(X>) (the Lindblad resonance) is canceled 
in the limit Sv^, — > 0. In general we will have to invoke 
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the assumptions (fl8|) and (|19p . though. Neglecting the 
small Sv,p term, the linearized continuity equation gives 



1 d . 

— — {roJ r ) = tAudp . 
r Or 



r5p = zil(^ 8E: 



Thus, 



Alu dr \m e j£l 
Integrating over dr and integrating by parts gives 

imQ/r qrpo 



rdrSp 



In the Lorentz gauge 



dr- 



(Alo) 2 m e -/Q 



(25) 



(26) 



(27) 



im 



8Ea, = 5$ + UijSAa 

r 

d 

8E r = — — 8$ + iuiSA r 
or 

SB Z = ™(r&4*) . 
r Or 



(28) 



For (Ar/r) 2 <C 1 we have the approximation 
r~ 1 d[r(d.../dr)] — d(...)/dr. Assuming the radial depen- 
dence (5$ oc exp(ifc r r) we obtain from the gauge condition 



ik r SA r H 6 A^ — icuSQ = 



(29) 



and 8A Z is negligible because of the symmetry of the 
problem. 8A r is non-zero since 8J r = poSv r + v r Sp 
and \5v^\ ~ 7~ 2 |Aw/ft||<5u r |. SAj can be computed from 
the Green's function (cf. appendix) to give 5A<p = v</,8$>. 
Fortunately, we only need 8E$. Note that 



SEtj, = -imr" 1 <5$/7 2 + O(Aw) 



IV. DISPERSION RELATION 



(30) 



The derivation of the Green's function can be found in 
the appendix 

/>oo 

5$(r) = 2ir 2 i r , dr , J m (uir < )H m {(jr > )Sp(r') (31) 
Jo 

The argument of the Bessel functions is assumed to be 
independent of r and r' in the important region 1 — 8 < 
r/ro < 1 + 8 which we will justify later. Thus, 



<5$(r) 



2tt ie m 



Hj 2 



' dr , n(r>)SS(r>) 
o r'(u> — mQ(r')) 2 



where Eq. [26] has been used. Since the right hand side 
of the last equation is independent of r $(r) has to be 
constant and we obtain with Q(r) rs 1/r 



2ir 2 ie 2 m 2 
Hj 2 



dr' 



n(f) 



r'(u> — m/r') 2 



(33) 



where the Bessel functions are evaluated at m(l — 
l/(27 2 )). The remaining integral can be evalu- 
ated if the Gaussian number density profile n(r) = 
no exp [— (5r 2 /2(Ar) 2 ] is replaced by a rectangle with 

width 2Ar and height noy/ir/2. 



dr'- 



n(r') 



r'(uj-m/r') 2 J ro{1 - S ) 
ln(wr' — m) 



r ^ 1+s) , , J^/2 
dr n — 



= n 



r'(u> — m/r') 2 

1 ro(l+<5) 



no 



7r 2r^8 



w 2 (wr' — to) 
to 



r (l-<5) 



2 a; (wro — m) 2 — w 2 Tq(5 2 



(34) 



The logarithm can be neglected because we ar e interested 
in the resonant case. With £ = 47re 2 no \Ar/ 2rg 8/m e j we 
obtain the dispersion relation 



1 = —TTt^Z"/ 



where Z = iJ m H m . Thus. 



1 



(Ad>) 2 - S 2 



Aw = ±^J8 2 - ir(Z/j 2 



(35) 



(36) 



where 



The Bessel functions can be expressed in terms of Airy 
functions for m 3> 1 

J m (z)= 2V3 m - 1/3 Ai(w;) (37) 
F ro (^) = - 2 J /3 m-VSgj^) ( 38 ) 

w = -2 1 /3 TO 2 / 3 (- - l) 

V TO / 

= (f) 2/3 [ 7 - 2 -2A^-2(£-l)] (39) 

with £ = r/ro which at the outer edge is £ — 1 ~ 8. 
Assuming m <C 27 3 = to 2 , 7~ 2 -C 1, |Aw| -C 7~ 2 
(assuming that the layer essentially exists in the region 

1 — 8 < r/ro < 1+5.) and 2 1//3 5 <C m~ 2 / 3 we can consider 

2 to be a pure function of m since |w| 2 <C 1. 



(0.347 + 0.200i)/m 2/3 



(40) 



Therefore, the Bessel functions do not depend on r and 
r' anymore. For \w\ 3> 1 the Airy functions can be ap- 
proximated to give 



1 



2ir\/w V rn 



2/3 



(41) 



but we additionally have to demand 8 <C 7~ 2 in order 
for Z to remain independent of r and r . For zero en- 
ergy spread the condition |Aw| <C 7~ 2 is equivalent to 
to <C 7 3 ^ 3 / 2 = mi. As we will discuss in Section V this 
condition may be too strict. 
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FIG. 2: Growth rates as a function of azimuthal mode m for 
£ = 0.02, 7 = 30 and various energy spreads. 



Some results are plotted in Fig. (2] For zero 
energy spread we recover the usual scaling relation 
3f(u;)/Q(ro) oc m 2 l z for m < 2 7 3 0, 0, d and 
3(w)/f2(r ) oc m 1 / 3 for m > 2 7 3 [| with the |io| > 1 
approximation of the Airy functions. The most dramatic 
consequence of a non-zero energy spread is the presence 
of a very sudden and steep cut-off. Pushing the cut- 
off to higher values of m requires increasingly small en- 
ergy spreads. For large energy spreads the growth rate 
3(w)/f7(Y(j) scales as m 1 / 3 if the \w\ <§; 1 approxima- 
tion is used. Retaining the Airy function the growth rate 
grows more slowly. According to Fig. [5] the scaling rela- 
tion for 5 = I/7 2 is in the order of to 1 / 4 before the cut-off 
is reached. This power law has been derived before in 
3 where the decoherence is due to betatron oscillations 
instead of a non-zero spread in the average angular veloc- 
ity. Note that these two results also agree with computer 
simulations that were carried out for Brillioun |§] flows. 



V. RESULTS AND DISCUSSION 

In order for the growth rate to vanish the expression 
inside the root must be real and non-negative. Using our 
approximations for the Bessel functions it is seen imme- 
diately that the former condition is satisfied if m > m,2- 
This explains why the drop off occurs at m — m,2 regard- 
less of the value of S as long as the remaining real part is 
positive. This is the case if the field reversal parameter 
does not exceed the critical threshold 



( > rwyS 



(42) 



If however this inequality is satisfied the expression inside 
the root becomes negative for m > mi and an unstable 
solution exists, but the sharp cut-off is replaced by Eq.l4*2l 
As discussed in the last section the azimuthal mode 
number m resulting in the largest growth rate is cither 
given by m = 2 7 3 or m = ^~ 1 S~ 2 whatever is greater. 
This guarantees that (A£>) 2 is either still complex or real 
and negative with both leading to an unstable mode. 
For our allowed parameter range 2 7 3 is typically larger. 
Thus, 



7- 1/2 - P 



(43) 



If the last inequality is not satisfied complex roots can 
only exist if Z is complex which is the case for m < no- 



where uj p is the nonrelativistic plasma frequency. 

Finally, let us consider small values of the azimuthal 
mode number m. So far we have have approximated 
the Bessel functions by Airy functions which makes them 
easier to compute especially for large orders. For small 
m this approximation cannot be justified and the Bessel 
functions have to be retained. In Figure [3] we solved the 
dispersion relation in the small m regime. Also, their 
arguments depend on Aw since m < mi. An accu- 
rate calculation of the growth rates of modes with az- 
imuthal mode nubers m is important for the determi- 
nation of the total power radiated, because the latter 
decreases with increasing m. Even for m = 1 Eq. [37] 
and [38] give excellent results, presumably because the 
argument of the Bessel functions is close to their or- 
der for 7 3> f and Z is independent of w if \w\ 2 = 

(f ) 2/3 [ 7 " 2 - 2Aw - 2(£ - 1)] < 1. The latter con- 
dition may be satisfied even if |A(D| violates our assump- 
tion \AQ\ < 7 ~ 2 . 

Looking at Eq. [34] we conclude that the quenching of 
the instability, i.e. the existence of the additional S 2 term 
in the denominator, is due to both the non-zero thickness 
of the layer and f2' ^ 0. However, it is important to note 
that the instability relies on the negative mass effect and 
would not exist in the absence of shear in our model. 
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APPENDIX: GREEN'S FUNCTION 

The Green's function for the potentials give 
(5$(r, t) = J dt'd 3 r' G(r-r',t- t') Sp(r', t') . 
6A(r, t)= J dt'd 3 r' G(r -r',t- t') <5J(r', t') , (A.l) 
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have for the electric potential, 
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FIG. 3: Growth rates as a function of azimuthal mode m for 
C, ~ 0.02, 7 = 30 and various energy spreads. The Bessel 
functions were retained in order to compute the growth rates 
in the low m regime correctly. 



where 



a 2 
at 2 



«3 



-4ir8(t)5(r) , G(k,w) 



47T 



k 2 



f OO p27T 

5$umk m (r) = 2 / r'dr' / ndn / da Sp^mk, (r') [ 



JO 



= 4-7T 

where 



/>oo />C 

/ r'dr' / 
Jo Jo 



kcIk 



K 2 — {uj 2 — k 2 ) 



(A.3) 



[■■] 



exp(zma) Jo{n[r 2 + (r') 2 — 2rr' cos a] 1 / 2 } 
k 2 — (uj 2 — k 2 ) 



where k 2 = k 2 + k 2 . Because to has a positive imaginary 
part, this solution corresponds to the retarded field. Also 
because Im(w) > 0, the n— integration can be done by a 
contour integration as discussed in [l6[ which gives 

/•OO 

^wmfc z (7') = 27r 2 i / r'dr'J m {kr < )H$(kr > )6p llli mkA' r ')i 
Jo 

(A.4) 

where k = (u> 2 — k 2 ) 1 / 2 , where r< (r>) is the lesser 

(greater) of (r, r'), and where Hm{x) — J m {x) +iY m (x) 
is the Hankel function of the first kind. Because 
(SAfaSJj,) = -(SA^SJx) cos (p + (5A y ,5Jy)sm<j>, we 
have instead of Eq. (|A.4|) , 



G(r, t) = 4tt J dio I d 6 k 



, 3 exp(ik • r — iut) 



k 2 



(A.2) 



where G is the Fourier transform of the Green's function. 
The "C" on the integral indicates an u>— integration par- 
allel to but above the real axis, Im(w) > 0, so as to give 
the retarded Green's function. 

Because of the assumed dependences of Eq. ([5]), we 



5^ mk *{r) = r8A" mK » = ir'ri / r'dr'6J" mK * (r') 

Jo 

■ J m+1 (fcr < ) J ff I ( ^ 1 (fcr > ) + J m _ 1 (kr < )H^_ 1 {kr > )(A.5) 

For to 3> 1, one can show that in this equation 

[Jm+iii^+i + Jm-i-ffm-i ]/ 2 = Jm Hm to a good ap- 
proximation. Equations (| A.4|) and (|A.5|1 are useful in 
subsequent calculations. 



[1] Manchester, R.N., & Taylor, J.H. 1977, Pulsars, (Free- 
man & Co.: San Francisco) 
[2] G. Benford and D. Tzach, Ap&SS 80, 307 (1981) 
[3] B. S. Schmekel, R. V. E. Lovelace, and I. M. Wasserman, 

Phys. Rev. E 71, 046502 (2005) 
[4] J. Arons and D. F. Smith, ApJ 229, 728-733 (1979) 
[5] R. C. Davidson, Physics of Nonneutral Plasmas, (Impe- 
rial College Press, London 2001) 
[6] B. S. Schmekel, Phys. Rev. E 72, 026410 (2005) 
[7] P. Goldreich and D. A. Keeley, ApJ 170, 463 (1971) 
[8] R. Buschauer and G. Benford, Mon. Not. R. astr. Soc. 

185, 493 (1978) 
[9] C. E. Nielson, A. M. Sessler and K. R. Symon, Proc. 
of International Conference on High Energy Accelerators 
and Instrumentation, pp. 239, (CERN, Geneva, 1959) 
[10] A. C. Entis and A. A. Garren and L. Smith, Proceedings 



of the 1971 Particle Accelerator Conference, pp. 1092- 
1096, (IEEE, Piscataway, New Jersey, 1971) 
[11] R. J. Briggs and V. K. Neil, Plasma Physics, 9, 209-227 
(1967) 

[12] A. A. Kolomenskii and A. N. Lebedev, Proc. of Interna- 
tional Conference on High Energy Accelerators and In- 
strumentation, pp. 115, (CERN, Geneva, 1959) 

[13] S. Heifets and G. Stupakov, Phys. Rev. ST Accel. Beams 
5, 054402 (2002) 

[14] O. Larroche and R. Pellat, Phys. Rev. Lett. 59, 1104 
(1987) 

[15] R. V. E. Lovelace and R. G. Hohlfeld, ApJ 221, 51 (1978) 
[16] G. N. Watson, A Treatise on the Theory of Bessel Func- 
tions, pp. 428-429, (Cambridge University Press, Cam- 
bridge, 1966) 



